Revisiting Absorbing Phase Transition in Energy Exchange Models 
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A recent study of conserved Manna model, with both discrete and continuous variable, indicates 
that absorbing phase transitions therein belong to the directed percolation (DP) universality class. 
In this context we revisit critical behaviour in energy exchange models with a threshold. Contrary 
to the previous claims [arAi?;:1102.1631, PRE 83, 061130 (2011)], our results indicate that both the 
maximal and minimal versions of this model belong to the DP class. 
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Absorbing state transitions (APT) [l|, Q in presence of 
a conserved field 0, 0] has drawn a lot of attention in 
the past decade. It is believed that critical behaviour 
of these systems is different from directed percolation 
(DP) the generic universality class of absorbing state 
phase transitions. The most important of these sys- 
tems is probably Manna model Q ; though originally in- 
troduced in context of self organized criticality Q, the 
conserved Manna model is well established in the litera- 
ture by now 0, H| . It is well known that the conserved 
Manna model shows an APT to infinitely many absorbing 
configurations when the density of the particles crosses 
a critical value <p c . There are several other model sys- 
tems which show similar APT in presence of a conserved 
field, namely conserved lattice gas conserved thresh- 
old transfer process Q etc. All these models show similar 
critical behaviour in 2D and higher dimensions; in fact 
it has been conjectured that absorbing phase transitions 
in stochastic models with infinite absorbing states and 
activity coupled to a non diffusive conserved field define 
a unique universality class 

Recently, however, it has been shown (gj that conserved 
Manna model actually belongs to the DP universality 
class. The apparent non-DP behaviour was due to the 
use of random initial conditions, which take unusually 
long time to reach the stationary state in these models. 
It was also shown that random initial condition shows 
an 'undershooting' in the decay of activity p a (t) which 
leads to erroneous determination of the critical point. 
This could be avoided by choosing suitable (natural) Q 
homogeneous initial conditions, which reach the steady 
state in reasonably short time allowing one to approach 
closer to the critical point. Corresponding critical expo- 
nents are found to be in good agreement with DP values. 
These results motivate us to revisit APT in some recently 
studied one dimensional systems where the dynamics is 
very similar to that of the Manna Model, which also ap- 
parently show non-DP critical behaviour 1^, 

The common features of these models U| are the 
presence of a conserved continuous local energy field, and 
an energy exchange dynamics; hereafter we refer to them 
as energy exchange models (EEM). The two versions of 
the EEM, maximal 11 1 and minimal [l(J models differ 
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FIG. 1. (Color online) Undershooting: Comparison of p a {t) 
versus t plot starting from random (green), deterministic (or- 
ange) and natural (blue) initial condition for both maximal 
(a) and minimal (b) models for a system of size L — 10 J .. The 
undershooting becomes more pronounced as one approaches 
the critical density e c . 

in the way the activity is defined locally. In the maxi- 
mal (minimal) version of the model a pair of neighbour- 
ing sites is called active when at least one of these sites 
has energy larger (smaller) than or equal to a predefined 
threshold value w. The active pairs reshare their total en- 
ergy randomly among themselves dynamically. Previous 
studies of EEM have claimed that these models undergo 
an absorbing phase transitions characterized by critical 
exponents different from those of DP. It is noteworthy 
that the dynamics of these models are similar to that of 
the continuous conserved Manna model (CCMM), where 
an active site relaxes by distributing all its energy ran- 
domly among the neighbours. Since CCMM is known 
to have DP critical behaviour Q, it is natural to expect 
that EEM models might also belong to DP class. In par- 
ticular we have observed that the activity density p a (t) 
in EEM models shows an undershooting when the ini- 
tial condition is random (see Fig. [1]). As it is known 
Q that undershooting in p a (t) leads to erroneous esti- 
mation of the critical point, we avoid these ill effects in 
the EEM models by using a regular homogeneous ini- 
tial condition. This leads to a more precise estimation 
of the critical point and critical exponents revealing that 
the true critical behaviour of APT in both maximal and 
minimal energy exchange models are in DP class. 

The Maximal Model: The maximal version of en- 
ergy exchange model is defined on a periodic one dimen- 
sional lattice with L sites labeled by i = 1, 2, . . . L; each 
site i containing a continuous energy Ei . A pair of neigh- 
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FIG. 2. (Color online) Critical point e c = 0.75243 for the 
maximal model is determined from p a {t)t a versus t. The inset 
shows log scale plot of p a {t) versus t; the slope of e = e c curve 
measured over last two decades corresponds to a = 0.1598(1). 

bouring sites is called active when at least one of the sites 
have energy larger than or equal to a predefined threshold 
w. Otherwise, i.e. when both these sites have energy less 
than w, the pair is called inactive. The system evolves 
following a random sequential update rule where a ran- 
domly selected active pair (i, i + 1) reshare their energies 



Ei^r(Ei + E i+1 ) 
E l+l (l-r)(-E( +E i+1 ). 



(1) 



Here r £ (0, 1) is a uniform random number. Clearly 
this dynamics is energy conserving, i.e. the average en- 
ergy e = Yli=i E% does not vary with time. Note, that 
similar conserving dynamics have been studied in differ- 
ent contexts of heat transport [l2[ and Econo-physics 13 1 
earlier. 

A configuration of the system where none of the pairs 
are active is called absorbing; it is easy to see that for 
e < w, the system has infinitely many absorbing configu- 
rations. However for e > w the system is active as there 
can not be any absorbing configuration for at least one of 
the sites must have Ei> w. This entails a possibility of a 
phase transition from absorbing to active phase when the 
average energy e is increased beyond some critical value 
e c , which is smaller than the threshold value w. Without 
any loss of generality we proceed with w = 1. Any other 
arbitrary w would only scale the critical point by a factor 
w, i.e. e c — > we c as the dynamics is unaltered under this 
scale transformation. 

The APT here is characterized by the average density 
of active pairs p a = (t^ ) , where Tj = 1 or depending on 
whether the pair (i, i + 1) is active or inactive. In the long 
time p a (t) saturates to a steady value, which is non zero 
only in the active phase and serves as the orderparameter 
of this transition. It has been claimed earlier that the 
critical behaviour of the model is different from DP. In 
the following we revisit the model carefully and show that 
the actual critical behaviour belong to DP. The primary 
reason for the apparent non-DP behaviour is possibly a 
small error in estimation of the critical point, which may 
be blamed to the use of random initial condition that 
shows undershooting. 

The usual random initial conditions do not contain 
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FIG. 3. (Color online) Estimation of /3 for the maximal model: 

(a) Plot of p a {i) versus t for various energy density e > e c . 

(b) Corresponding saturation values plotted against e — e c in 
double logarithmic scale; the slope corresponds to /3 = 0.283. 

the natural correlations existing in the stationary state. 
In particular the density profile in a random configu- 
ration is disordered as opposed to the stationary pro- 
file which is essentially flattened by the underlying dif- 
fusive dynamics. As discussed in Ref. Q, one should 
use the reactivated stationary configurations, called nat- 
ural initial conditions for studying decay of activity in 
models with additional conserved field. However, gen- 
crating natural initial condition takes substantial com- 
putational time which can be avoided in these continu- 
ous models by artificially creating a homogeneous density 
profile {Ei} = {2e, 0, 2e, . . . }. This deterministic initial 
configuration is fully active with p a = 1 when e > 1/2. 
Fig. Ufa) shows a comparison of the decay of p a {t) in the 
super critical regime starting from three different initial 
conditions : random (green), deterministic (orange) and 
natural (blue). Both the natural and deterministic initial 
condition do not show any undershooting as opposed to 
the random initial condition. In this work, we choose to 
study APT using deterministic initial condition as it is 
computationally inexpensive. 

Critical point and a : At the critical point e c the ac- 
tivity decays as a power law 



Pait) ~ t" 



(2) 



One can estimate the critical point e c and exponent a by 
plotting p a (t) versus t for various values of e and looking 
for a power law decay. This estimate can be verified from 
the plot of p a (t)t a against t; the curve corresponding to 
e = e c would remain constant in the long time limit. This 
procedure is illustrated in Fig. [2] which gives an estimate 
e c = 0.75243(3). The log scale plot of p a (t) at e c gives an 
accurate estimate of the critical exponent 

a = 0.1598(1). 

This value 1 of o is a bit lower than the previous estimate 
0.19 [ll| and is in good agreement with the DP value 
(see Table I). The difference in a here originates from 
the current estimate of e c which is slightly higher than 
the previously reported [llj value 0.7508. Note that the 
system size used for the current work (L = 10 5 ) is much 
larger. This correction of e c is also expected to modify all 
the other critical exponents of the system substantially. 
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FIG. 4. (Color online) Finite size scaling for the maximal 
model : p a (t)L l3l ' ± for different L plotted against t/L z col- 
lapse for z = 1.58 and fi/v±_ — 0.26. The inset shows the 
corresponding unsealed data. 

Stationary state and (3 : In the active phase the ac- 
tivity saturates to some finite value p a which vanishes 
algebraically as one approaches the critical point, 



p a ~ (e - e c ) 



(3) 



Here /3 is the order parameter exponent. Figure OUJa) 
shows p a {t) as a function of i for various values of A = e— 
e c in the supercritical regime. Corresponding saturation 
values are plotted against A in log-log scale (see Fig. 
[3ft)) ) ; the slope of the resulting straight line gives us an 
estimate 

/3 = 0.283(4). 

Again this value of (3 is much lower than the previous 
estimate and is consistent with the DP value Pdp = 
0.2764. Note that the present estimate is done for very 
small e — e c < 10~ 3 which was not achieved earlier. 

Finite size scaling and z: Now we turn our attention 
to the finite size scaling. For a finite system, the decay of 
p a {t) at the critical point is expected to follow a scaling 
form 



Pa (t) = L-^Q(t/L z ), 



(4) 



where z is the dynamic exponent. Thus, one expects that 
at the critical point e = e c , activity density p a (t)L^I Vl - 
for different values of L would collapse on to a single 
function when plotted against t/L z . This data collapse, 
shown in Fig. 0] for systems of size L = 2 8 , 2 9 . . . 2 12 , 
yields an estimate 



2 = 1.58(2) ; - =0.26(1 



V _L 



which are, again, in excellent agreement with the corre- 
sponding DP exponents. 

In summary, all the critical exponents a, ft, z and j3/v± 
for the maximal model agree quite well with the corre- 
sponding DP values (see Table U). On the basis of these 
estimates we conclude that, contrary to the earlier claim, 
the maximal energy exchange model indeed belongs to 
the directed percolation class. The errors in the previ- 
ously reported [111 ] exponents were caused by the inaccu- 
rate estimation of the critical point resulted from the use 




FIG. 5. (Color online) Critical point e c for the minimal model 
: Pa(t)t a versus t curve for a system of size L = 10 5 remains 
constant in the long time regime for e c = 1.24096. The inset 
shows estimation of a = 0.1595 from p a {t) vs. t plot. 

of small system size and random initial condition (which 
suffers from undershooting). 

The Minimal Model : In this section we revisit the 
minimal version of the model studied recently jioj ]. We 
restrict ourselves to one dimensional version of the model 
only. No particular difficulty, except the requirement of 
large computational time, is expected in studying the 
model in higher dimensions. Again the minimal model 
is known to have same mean field exponents as those of 
DP, the departure from DP, if at all there, should be 
prominent in one dimension (ID). 

In the ID minimal model, each site i has energy Ei, 
which is a continuous variable. A neighbouring pair of 
sites + 1) is defined to be active when at least one of 
the sites i or i+1 has energy less than or equal to a thresh- 
old value w. The active pairs in the model evolve accord- 
ing to the energy conserving dynamics |T]). It is clear 
that any configuration having average energy e = J] Ei 
smaller than w is active. Configurations where all the 
sites have energy larger than w are absorbing; the sys- 
tem is likely to be trapped in an absorbing configuration 
for e ^> w. Thus, a transition from an active phase to 
an absorbing state is expected as the average energy is 
increased beyond some critical value e c > w. 

We must mention that the earlier study of the mini- 
mal model [l(J was done at a fixed e = 1 and varying the 
threshold w. This resulted in a critical point w c = 0.810; 
in other words in the (e, w) phase plane (1, w c ) is a crit- 
ical point. It is evident that (A, Xw c ) is also a critical 
point for any arbitrary A as both the conditional state- 
ment Ei < w and the dynamics ([1]) are invariant under 
a scale transformation (Ei — >• XEi,w — > Xw). This im- 
plies that the critical line in the (e, w) phase plane is 
a straight line with unit slope passing through (l,w c ). 
Clearly (l/w c ,l) is also a critical point and it can be 
reached keeping w fixed; we choose to study the critical 
behaviour of the minimal model along the line w = 1 by 
varying e, expecting the transition at e c = l/w c . 

We proceed further by defining the density of active 
pairs p a as the order parameter. At the critical point 
Pa{t) ~ t~ a . The decay of p a {t) is studied starting 
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FIG. 6. (Color online) Estimation of j3 for the minimal model : 
(a) p a (t) versus t for different e in the supercritical regime, (b) 
Corresponding saturation values as a function of A = e c — e; 
the slope of the solid line corresponds to j3 = 0.2749. 

from the spatially ordered homogeneous initial condition 
{Ei} = {2e, 0, 2e, 0...}. Unlike random initial con- 
dition this homogeneous initial configuration does not 
suffer from undershooting (Fig. [lib)). 

Repeating the standard procedure for the determina- 
tion of the critical point, discussed in section II, we plot 
Pa(t)t a versus t for different e in Fig. O This curve 
becomes flat at large t for 

e c = 1.24096(3). 

The estimate of critical point is slightly higher than the 
previously reported value l/w c — 1.23456 [l(J ■ Con- 
sequently, a careful study of the critical exponents in 
large system size shows a considerable departure from 
earlier estimates and agree remarkably with DP univer- 
sality class. 

The slope of the curve p a (t) versus t in log scale (over 
last two decades) at the critical point e c , shown in Fig. 
gives an estimate of the exponent 

a = 0.1595(4), 

which is in strikingly good agreement with aup- 

Next, we calculate the order parameter exponent f3 
from the steady state density of activity p a . Figure Eta) 
shows the saturation of p a (t) for various values of e rang- 
ing from 1.24080 to 1.21540. The stationary densities p a 
are plotted against e — e c in log scale in FigJHJb) to obtain 
an estimate of (3; a linear fit for last two decades gives 

P = 0.2749(8), 

which, again, matches well with ajjp. 

The dynamical exponent z and /3/v± can be obtained 
from the finite size scaling, as at the critical point the 
activity density of finite systems scales as Eq. (|4]). In 
Fig. [7] we have shown the decay of p a (t) (inset) for 
L = 2 s , . . . , 2 12 . These curves have been collapsed onto a 
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TABLE I. Comparison of critical exponents. 




FIG. 7. Finite size scaling for the minimal model according 
to eq. (|4]). The best collapse corresponds to z = 1.54 and 
P/vx = 0.25. The inset shows the unsealed data. 

single scaling function in Fig. [7] using 

z= 1.54(2) ; ^ =0.25(1). 

Both the estimates are in good agreement with the cor- 
responding DP values. 

To summarize, the critical exponents of the minimal 
model also imply that, contrary to previous claims fioj ]. 
this APT belongs to DP. We should mention that these 
differences in exponents arc not in an ywa y related to the 
fact that the orderparameter in Ref. [10j was chosen dif- 
ferently. Instead of density of active pairs p a , the average 
density of sites having Ei < w was used there. This or- 
der parameter p = (s,), where Sj = 1 (sj = 0) if Ej < w 
(Ei > w), is related to p a as 

Pa = P+ (Sj(l - S i+ l)). (5) 

Since p a is nonzero only when p ^ and they are dimen- 
sionally identical, the critical point and corresponding 
exponents would be same in both cases. 

In conclusion, we have revisited absorbing phase tran- 
sition in the energy exchange models in one spatial di- 
mension. With extensive Monte-Carlo simulation on 
large systems, along with homogeneous initial conditions, 
we establish that both the maximal and the minimal ver- 
sion of the model belong to DP universality class; the 
critical exponents are summarized in Table [I] 

With this study we want to emphasize the fact that one 
has to be very careful in exploring absorbing phase tran- 
sition in presence of a conserved field. The conserved 
field imposes additional constraints on the evolution of 
the system and introduces a long time scale. Random 
initial conditions may need very long time to saturate, 
producing unwanted transient effects like undershooting. 
These undesirable features may lead to erroneous esti- 
mation of critical point and hence inaccurate critical ex- 
ponents. One must take care of these ill effects while 
studying APT in a system with a conserved field. 
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